Aging in Ki_ x Li z Ta03: a domain growth interpretation 
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The aging behaviour of the a.c. susceptibility of randomly substituted Ki-^Li^TaOs crystals 
reveals marked differences with spin-glasses in that cooling rate effects are very important. The 
response to temperature steps (including temperature cycles) was carefully studied. A model based 
on thermally activated domain growth accounts for all the experimental results, provided one allows 
for a large distribution of pinning energies, in such a way that 'slow' and 'fast' domains coexist. 
Interesting similarities with deeply supercooled liquids are underlined. 
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The a.c. susceptibility of 'glassy' materials becomes 
time dependent at low temperatures: this is the aging 
phenomenon [nj. Such a behaviour has been observed on 
spin-glasses , on disordered dielectrics such as ammo- 
nium pcrchlorate 01 or Ki-^Li^TaOa (klt) (which is an 
orientational glass [||]6|), and more recently on a deeply 
supercooled liquid (glycerol) 0. All these systems share 
common properties, but also important differences. The 
prominent features are the following: 

• When the system is cooled to a certain temperature 
T\ below the glass transition Tq, its frequency dependent 
response (which in the present paper is the real part of 
the dielectric constant e') depends on the time At since 
the quench. This dependence can be parametrized as pi : 



e'(u,t) = e' sr (u) + f(At)g(u) 



(1) 



where /(At) is a (slowly) decreasing function describing 
the aging part and e' ST (uj) is the stationary part of the 
response. In spin glasses, the functions f(u) and g(u) 
both behave similarly (as u~ b j^Jj)), which means that 
the aging part of the response obeys approximatively a 
simple u) ■ At scaling (but see j||l|]). On klt crystals, 
g{u>) is nearly independent of frequency, at least for tem- 
peratures not too close to T , while /(At) behaves as a 
power law with a small exponent j^] or a logarithm, see 
Fig. 1. The situation is intermediate for glycerol, where 
the variation of g(u) is somewhat weaker than that of 

/(«) §■ 

• The stationary part of the susceptibility is nearly in- 
dependent of the cooling rate TZ in spin glasses || . This 
is in striking contrast with the case of klt, where e' ST (oj) 
markedly depends on the cooling rate (see Fig. 1), de- 
creasing roughly as log TZ, as the cooling rate is reduced. 
In other words, the apparent asymptotic value of the di- 
electric constant depends on the history, and thus cer- 
tainly cannot be associated to an 'equilibrium' response 
function. Several cooling histories were probed in ref . || , 
with the conclusion that it is essentially the time spent 
around To which determines the value of e' ST (u>). 
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FIG. 1. Plot of e'(u/2n = 100 KHz, t) for klt 
(x = 0.025) as a function of the time At = t — t\ 
spent at temperature T\ = 4.2K, for various cooling rates: 
U = 0.045, 0.024, 0.011, 0.0058 K/s. The slowest cooling 
rate corresponds to the lowest value of the apparent asymp- 
tote e' ST (u>). Insets: EXP: same data, in a log-linear scale; 
rfim: numerical simulation of the time evolution of the do- 
main wall density in the 3-D Random Field Ising Model, for 
different numerical cooling rates Q . 



• Upon a negative temperature cycling T\ — » T2 = 
T\ — AT — » Ti, the a.c. response function shows both 
rejuvenation and memory, in the sense that new out of 
equilibrium processes are induced during the low tem- 
perature interval, but are erased once the temperature 
is back at T%: dynamics resume as if the intermediate 
regime had not taken place. Memory is nearly perfect in 
spin-glasses when AT > IK, while in klt (or glycerol) 
there appears an interesting 'overshoot' - see Fig. 2 be- 
low - before the dynamics corresponding to the initial 
temperature can resume. The amplitude of this over- 
shoot grows with the time ti spent at temperature T 2 . 

The aim of this paper is to show that a picture based 
on domain growth, where domain walls are pinned by 



1 



static impurities, does account very well for all the exper- 
imental data on klt, provided one takes into account the 
fact that the growth rate is actually very inhomogeneous, 
which leads to the existence of both 'slow' and 'fast' do- 
main walls. The relation with the 'droplet picture' for 
spin-glasses will also be discussed in detail. Interestingly, 
the similarities between klt and glycerol suggest that 
some sort of pinned domain growth, with a character- 
istic length scale, might also be present in supercooled 
liquids. 

When lithium is absent, KTa03 is known to be an 
'incipient' ferroelectric at T = 0, where quantum fluc- 
tuations prevent total ordering. Thus, the ferroelectric 
correlation length £(T) is expected to be large at low 
temperatures: assuming that £ a one finds that 
f(T) cx 1/VT. When lithium ions are randomly substi- 
tuted, they create large local dipolcs which are known to 
freeze individually around T = AOK, due to the energy 
barrier for dipole reorientation. Below this temperature, 
the lithium ions can thus be seen as giving rise to a ran- 
dom static electric field, which act as a pinning field for 
the ferroelectric domain walls. The growth of the ferro- 
electric order towards its equilibrium value £(T) is thus 
strongly impeded. The model we have in mind in thus 
the random field Ising model (rfim) above its ferromag- 
netic transition, but close to it so that the equilibrium 
domain size is somewhat larger than the lattice spacing 
a. A numerical simulation of the three dimensional rfim 
with different cooling rates actually reveals qualitatively 
similar features (Fig 1 - inset) ||. 

In the presence of pinning, several arguments [ ^0|JTT| 
suggest that the typical energy barrier which impedes 
upon the motion of a domain of size R is of the order 
of E{R) ~ T (~ ) where 9 ~ | for nearest neighbour 
interactions in three dimensions and T is an energy scale 
which depends both on the random fields and on the sur- 
face energy of the domain walls. Note that dipolar fields 
are long range, so that 9 might have a somewhat larger 
value (9 ~ 2). The effect of pinning becomes noticeable 
when E(R) ~ kT. For large times, the growth law is 

thus governed by thermal activation t ~ To exp E ^ , 
where tq is a microcopic time scale, and the Boltzmann 
constant k has been set to one. One thus obtains |]l0| , |lT| : 

R(t) ~ (^ log(t + e r <^ 6 R = R(t = 0) (2) 

where from now on we measure t in units of To and R 
in units of a. Eq. (^J) holds until the equilibrium size 
£(T) is reached, beyond which domains no longer grow. 
Taking the derivative of Eq. (0) , one obtains the growth 
rate T(t, T), in which the temperature can now be consid- 
ered as time dependent, for example as T(t) = Tq — IZt. 
Noting that at t = domains are in equilibrium at Tq 
(i.e. Rq — £(2q)), one can integrate back T(t, T(t)), and 



find the domain size as a function of time at tempera- 
ture T\. In order to compare with experimental results, 
we make the natural assumption that the domain walls 
contribute to the dielectric susceptibility proportionally 
to their total surface. Per unit volume, this leads to an 
excess susceptibility given by: 

^'{u,t)<x{R(tY)/{R{tf)g(u) (3) 

where g{oj) is the (frequency dependent) wall mobility 
and (...) denotes an average over all the domain sizes 
present in the sample. 

Numerically, however, this simple model fails to repro- 
duce the experimental results in the following sense: T 
is either too small, and one observes a decay of e', which 
can be slow but towards an asymptotic value which is 
independent of the cooling rate, or too large, in which 
case e' indeed strongly depends on the cooling rate but 
does not relax at all at T\ (at least over the experimental 
time window) . The simultaneous observation of the two 
effects actually requires some fast growing domains coex- 
isting with slower ones, corresponding to different local 
barrier heights T. This is actually expected: since the 
pinning field is random, Eq. (Q) can only describe a typ- 
ical behaviour, with large fluctuations corresponding to 
particularly efficient (or inefficient) pinning regions. We 
shall thus assume that T has a rather broad distribu- 
tion. Let us first consider the situation at time t\ when 
the temperature T\ is first reached (i.e. To — lZt\ =T\). 
There, three types of domains can be distinguished: 

• 'Fast' domains, which have been able to remain in 
equilibrium during the whole cooling phase, and follow 
the evolution of the equilibrium length £. Their size is 
thus £(Ti) and they no further evolve after t\. Imposing 
that T(t, T) > d£/dt, one can check that fast domains 
correspond to T p < £T 1 1+e/2 ' /T° /2 ', where T p = T£(T ) e 
is the 'initial' pinning energy and £ is approximately given 
by HTi/KTo}. 

• 'Slow' domains, which are in equilibrium in the first 
stages of the cooling process, but fall out of equilibrium 
before the temperature reaches T±. This occurs when 
the domain growth rate T(t, T) becomes less than d£/dt, 

which occurs at a temperature (T 2 Tq /I 2 } 1 . This 

corresponds to £T^ +e/2 /T® /2 <T p < £T . 

• 'Frozen' domains, which fall out of equilibrium as 
soon as T < To because of a large local pinning energy 
T p > £Tq. Of course, these domains are not really frozen 
but only evolve on extremely long time scales. 

Collecting the contribution of these different domains 
to (R(ti) 2 ) / (R(ti) 3 ) , and assuming that the scale T* of 
the distribution of pinning energies T p is much larger 
than £Tq, we find |L9| that the contribution of the 
walls to the dielectric susceptibility is proportional to 
1 — crfXo/T*, where a is a certain function of Ti/Tq 
p9[ . Hence, this model predicts that the excess dielec- 
tric susceptibility is linear in In 1Z, as indeed reported 
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in pj. The model also allows us to compute the subse- 
quent evolution of e. Physically, after a time t± + At, 
the 'fastest' slow domains have managed to reach £(Ti), 
while the others slowly grow. Adding the different con- 
tributions, we find that, for large At (more precisely for 
£{T 1 /T )( 1 ~ 9 ')/ 2 < In At < T*/T ), the dielectric con- 
stant decreases as — j3T$/T* In At, where is another 
function of Ti/T (but independent of 1Z). This loga- 
rithmic decay is again in good agreement with the exper- 
imental data (at least for large enough At): see Fig 1. 
Of course, for exponentially long times, all domains (even 
the 'frozen' ones) reach their equilibrium size the 
point is that the logarithm is such a slowly varying func- 
tion that, even on the rather long time scales investigated 
here (3 10 5 seconds), the dielectric susceptibility appears 
to asymptote a cooling rate dependent value e' ST (iu). 
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FIG. 2. Effect of a negative temperature cycle on 
e'(oj/2-K = 100 KHz, t) for klt (a; = 0.011). The first pe- 
riod is spent at temperature T\ = 18.4 K during 40000 s., 
then the sample is cooled to 4.2 K for a time t-z = 40000 s., 
and heated back again to T\. The 'overshoot' can be seen 
at the beginning of the third period, where e'(u>,i) has a non 
monotonic behaviour. 



be. However, since the size of these domains is large, the 
barriers they have to jump are now high at T\ and their 
shrinking is very slow. A more efficient process then oc- 
curs: internal nucleation of smaller domains. These nu- 
clei grow according to Eq. (|^) from an initial size of 
the order of the lattice parameter a, until they reach the 
equilibrium domain size £(Ti). This growth, rapid at the 
beginning, generates new walls which induce an increase 
of the dielectric constant. This accounts for the 'over- 
shoot' observed just after the temperature step. Second, 
the 'slow' domains which have been nearly unaffected by 
the temperature decrease, resume their dynamics in per- 
fect continuity with the first period of time. Concerning 
the 'overshoot', its amplitude A(t 2 ) depends on the time 
t 2 spent at low temperature, since the number of net 
nucleation of domains is proportional to the volume oc- 
cupied by domains which have grown beyond £(Ti). We 
thus expect Afo) to grow slowly with time (since R(t 2 ) 
only grows logarithmically) and then to saturate for very 
long times. This is indeed what one observes experimen- 
tally. The overshoot contribution can be parametrized 
as A{t2) exp — {At/ri), where A(t 2 ) grows slowly with 
t 2 , and decreases with T±, a feature which our model also 
predicts []l9| . T\ , on the other hand, only weakly depends 
on t 2 or T\ (see Figure 3). 
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The interpretation of the cycling experiments, shown 
in Fig 2, follows very similar lines. The main ingredi- 
ent is again the existence of the 'fast' and 'slow' domains 
discussed above. When the temperature is reduced from 
T\ to T2, the equilibrium length £ suddenly increases. 
Hence, the 'fast' domains are driven out of equilibrium 
and must restart growing. Since the temperature has de- 
creased, this is now a slow process, which accounts for the 
fact that e' decreases much like after the initial quench. 
(Note that the large positive instantaneous change of e' 
is due to the bulk contribution which is itself tempera- 
ture dependent). The 'slow' domains, which were already 
slow at T\, are completely frozen at T 2 and hardly move. 
Now, when the system is heated back to T x , two things 
happen. First, the 'fast' domains now have to shrink 
back to their equilibrium size. They have grown beyond 
£(Ti) during their stay at lower temperature and conse- 
quently, the domain wall density is smaller than it should 
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FIG. 3. Amplitude of the overshoot A(t2) as a function of 
the time ti spent at T2 = 4.2K, in log-linear scales, and for 
different temperatures T\. Note that A increases roughly as 
In £2 and decreases when Ti increases. Inset: Dependence of 
the relaxation time of the overshoot ri on t^- The j/-axis scale 
is - 2000 seconds. 



Finally, we would like to discuss the differences be- 
tween domain growth in a random field like system, 
which we argue to be a good model for klt crystals, 
and 'droplet' growth in spin-glasses |l2|-|l4j], which has 
been advocated to be the relevant picture (as opposed 
to models inspired from Parisi's mean field 'hierarchical' 
solution 1 15 16|,|l8 ], or exactly soluble dynamical mean- 
field models Jl7|,[l|). First of all, domain growth cannot 
be trivial in spin glasses, otherwise strong cooling rate 
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dependence would be seen, for example in the apparent 
asymptote of x'(w,i), as is observed in KLT (see Fig. 1). 
This is avoided by arguing that the growing 'phase' de- 
pends chaotically on temperature ]T^-|T4]|, i.e. the struc- 
ture of the equilibrium phase changes altogether when 
temperature is varied. Hence, the growth of the domains 
at temperature T + AT is useless to bring the system 
closer to equilibrium at T. The rejuvenation effect seen 
in klt was attributed above to a change of the finite 
equilibrium correlation length £ when the temperature is 
reduced, which leads to an overshoot effect not observed 
in spin-glasses. Conversely, in the droplet picture, the 
equilibrium size of the domains is infinite ||l3|] . The inter- 
pretation of rejuvenation is then related to the fact that 
domains of the new phase (i.e. the one stable at temper- 
ature T2) do grow at the expense of the old phase (the 
one stable at temperature Ti). However, this must be 
also compatible with the observed perfect memory upon 
reheating. How this is possible is still very much a sub- 
ject of debate (see e.g. 16§), but the conclusion that 
different phases must somehow coexist (at least for finite 
times) appears difficult to avoid. 

To summarize, we have argued that the observed ag- 
ing effects and strong history dependence of the dielec- 
tric susceptibility in klt crystals could be understood 
in terms of slow, inhomogeneous ferroelectric domain 
growth in the presence of random pinning fields. An 
important aspect, needed to interpret the 'overshoot' ob- 
served after a positive temperature jump, is that the equi- 
librium correlation length £ is finite. It would be inter- 
esting to compare our results with similar experiments 
on well characterized random field systems, in particu- 
lar in the ordered phase where £ = 00. We have dis- 
cussed the similarities and differences with other systems, 
such as spin-glasses or supercooled liquids. In particular, 
no overshoot has been reported in spin glasses for large 
enough temperature jumps Ti — T\, while the data on 
glycerol is quite similar to the one discussed here, and 
raises the interesting possibility that some kind of do- 
main growth, with a finite equilibrium size, might also 
be relevant in supercooled liquids. 

We thank J. Kurchan for having suggested the analogy 
with random field systems and M. Mezard, R. Leheny, P. 
Nordblad and E. Vincent for many useful discussions. 
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